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Abstract 

As a step toward satisfactory understanding of the quantum dynamics of Dirichlet 
(D-) particles, the amplitude for the basic process describing the scattering of two quan- 
tized D-particles is computed in bosonic string theory. The calucluation is performed and 
cross-checked using three different methods, namely, (i) path integral, (ii) boundary state, 
and (iii) open-channel operator formalism. The analysis is exact in a' and includes the 
first order correction in the expansion with respect to the acceleration of the D-particles. 
The resultant Lorentz-invariant amplitude is capable of describing general non-forward 
scattering with recoil effects fully taken into account and it reproduces the known result 
for the special case of forward scattering in the limit of infinitely large D-particle mass. 
The expected form of the amplitude for the supersymmetric case is also briefly discussed. 
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1 Introduction 



A recent proposal on microscopic formulation of M theory [|l], [| suggests that the Dirichlet 
(D-)particles, which are thought to describe solitonic collective excitations of string 
may actually play a much more fundamental role. In that approach, at least in the 
11 dimensional infinite momentum frame, D-particles are the basic degrees of freedom 
and extended objects such as membranes, strings, etc., arise as collective excitations. A 
certain amount of evidence exists in support of this aspiring conjecture, essentially in 
low energy domain. It is extremely important to check if the details work out, especially 
whether there would be higher derivative corrections to the proposed action. To answer 
this question, obviously we need to know much more about the quantum interaction of 
D-particles. 

The conjecture above was in part motivated by the effective low energy description of 
the D-brane interactions in the form of super- Yang-Mills theory on the worldvolume 
Although conceptually independent of M-theory, some hints about the hidden 11 dimen- 
sional feature have been uncovered 000, suggesting that extension of this approach is 
a viable road to understanding of the underlying fundamental theory. Another intriguing 
feature that appears in this approach is the apparent non-commutative nature of the co- 
ordinates describing the D-branes. Although natural from the gauge theory point of view, 
its physical meaning and how fundamental it actually is yet to be clarified. Again more 
detailed understanding of the quantum dynamics of D-branes is needed to make further 
progress. 

With this background, we present in this article a computation of the amplitude for 

scattering of two quantized D-particles with finite mass in bosonic string theory. The 

result, valid when the accelerations of the D-particles are small in their mass scale, is 

nevertheless exact in a' expansion and describes general non-forward scattering with recoil 

effects fully included. In the special case of forward scattering in the limit of infinitely 

heavy classical D-particles, the amplitude correctly reduces to the known expression ||. 

In fact, the structure of our result is formally quite similar to the one for the special case 

so that there exists a simple rule to go backwards, i.e. to go from this special case to our 

general case of quantum D-particles. Applying this rule, we will write down the result 

expected in the superstring caseQas well in the final section. Our result should serve as 

a firm knowledge to be checked against in any proposal for the fundamental microscopic 

theory containing D-particles. 

^dn superstring case, satisfactory analysis should include the processes involving fermionic superpart- 
ners of (bosonic) D-particles and this is yet to be performed. 
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The actual calculations are perfomed and cross-checked by three different methods, 
namely (i) the path integral, (ii) the boundary state, and (iii) the operator method in 
the open string channel. The emphasis is on the path integral method, which is concep- 
tually most complete in formulating the problem and hence is capable of computing the 
corrections in powers of acceleration systematically. It is a non-trivial extension of the 
formalism that we developed recently0 for the scattering of closed string states from a 
quantized D-particle. We will see that all the details work out consistently , including 
the first order correction^ that we will take into account in this article. As it will become 
clear later, the other two methods, although not suitable for computing corrections in 
powers of acceleration, have their own advantages. The boundary state method turned 
out to be most efficient for the purpose of computing the lowest order amplitude. On the 
other hand, the operator method reveals a deformation of the spectrum of an unusual 
type in the open string channel and the transitions among these excitations. Deeper un- 
derstanding of this phenomena may be quite important in connection with the effective 
gauge theory. 

In order to explain three different methods, each of which contains subtleties both 
technical and conceptual, in a reasonably self-contained manner, this article has become 
somewhat long. However, we have organized it in such a way that a hurried reader, if so 
desires, may consult section 2 on the path integral method, which by itself is complete, 
and skip sections 3 and 4 describing the other two methods. Here is the content of the 
rest of the article to follow: 

Section 2. is devoted to the path integral method. In subsection 2.1, after describing 
how DO-DO system should be characterized, we carefully analyze the boundary conditions 
when the D-particles follow arbitrary trajectories. They are then expressed in a suitable 
orthonormal coordinate system, to be used throughout the paper. Path integral over 
the string coordinate itself is performed in subsection 2.2. The complicated boundary 
conditions are implemented by devising a useful trick, and the lowest order amplitude as 
well as the Green's functions needed for the calculation of the corrections are computed. 
In subsection 2.3 the corrections, to first order in the acceleration, which arise from the 
extendedness of the boundaries are computed. The divergence produced in this process is 
consistently absorbed by the renormalization of the trajectories. Finally, in subsection 2.4 
we quantize the D-particle trajectories themselves and obtain the desired quantum am- 
plitude. It is then compared with the known result in the special limit already mentioned 

2 In this article, "order" invariably refers to the one in the power expansion with respect to the 
acceleration of the D-particles. 
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above. 

In section 3, we briefly describe how the lowest order amplitude is reproduced using 
the boundary state representation of the interaction vertex developed previously [p~Q|| P| . 
After writing down the appropriate vertex states in 3.1, we sketch the actual computation 
in subsection 3.2. A powerful normal-ordering formula is developed (a proof is sketched 
in the Appendix) and used to calculate the Green's functions as well as the amplitude. 

The computations performed in section 2 and 3 are essentially from the closed string 
channel. In section 4, we develop the operator formalism in the complimentary open string 
channel. First in 4.1, we make the modular transformation of the amplitude previously 
computed and obtain the form to be reproduced in the open channel. Upon performing 
the operator quantization in 4.2, we will find that the structure of the Hilbert space 
including the energy spectrum is rather unusual. We will give a careful analysis of this 
structure, give physical interpretation, and compute the relevant trace in a reliable way. 
The end result is exactly the one expected. 

Finally, in section 5, we make several important remarks that emerge from our work. 
One among them is about the expected form of the amplitude in superstring theory. We 
will spell out how, by applying the aforementioned rule, the result of || for the special case 
can be promoted to the amplitude for quantized D-particles. Another remark concerns 
the existence of infinite number of diagrams contributing to the scattering in the high 
energy regime even for small string coupling. 

2 Path Integral Approach 

2.1 Characterization of DO-DO system 
The setup 

A D-particle (DO-brane) is a point-like object which can emit and absorb a closed string. 
It does so in such a way that, when looked at from the (dual) open string channel, the ends 
of the open string lie somewhere on its worldline. Thus the simplest diagram describing 
two such D-particles interacting with each other looks like the one dipicted in Fig.l. 
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Fig.l Basic process for DO-DO scattering. 

In the open string language, the relevant topology of the string worldsheet is an annulus 
E, which we take to be one with the inner radius R\ and the outer radius i?2- The 
boundary circle at each end will be parametrized by the angle 9. We take the inner 
boundary to be mapped onto the worldline of a D-particle, parametrized by while 
the outer boundary is embedded into the worldline of the other D-particle, described by 
fliitz)- Until we come to the subsection 2.4, where we quantize the D-particles, we take 
fiifi), i = 1, 2 to be arbitrary yet fixed time-like trajectories^ . As the ends of the open 
string may terminate anywhere on the worldlines, the Lorentz covariant constraints at 
the boundaries should be of the form [11 

X»(r = R i ,9) = ft(ti(9)), (2.1) 

where X M (r, 9) denote the open string variables in the worldsheet polar coordinates and 
ti(9) are arbitrary functions describing the embedding. This means that in the path 
integral formulation we develop in this section, we must integrate over X M (r, 9) in the 
bulk and over ti(9) at the boundaries. Thus the relevant amplitude is 

v(/i,/ 2 ) = / p^(r, n^(^) n ^(^(^, ^) - /r^w)e- sw , (2.2) 

i i 

where 

S[X] = J^dPzd a Xi*d a Xp (2.3) 

is the open string action^. 
Scheme for i; (#)-integration 

To give a feeling for how we will compute the amplitude above, we need to briefly recall 
the scheme for ^(^-integration explained in || (for the single D-particle case). We first 
3 There will be a slight condition required on /f (U) for consistency later. 

4 We use Euclidean worldsheet and Minkowski target space with the metric r\^ LV = diag (— , +,+,••■,+). 
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split the integral over ti(0) into the one over the ^-independent mode, denoted by ij, 
and the rest over the non-constant modes, and then, to retain general covariance, expand 
//* in terms of the geodesic normal coordinate Q(9) around ft(U)- It can be written 

in the form 

ft(U(0)) = ft(ti) + fmCi(0) + ^(9) : (2.4) 
where Q?(e) = \k^? 

Here, a dot stands for a t-derivative and hi(ti), Kf{ti) and P^ v {ti) are, respectively, the 
one-dimensional induced metric on the trajectory, the extrinsic curvature and a projection 
operator normal to the trajectory. They are given by 

hi = f?U, (2.6) 

K? = ft-~ft = PrU, (2-7) 
= rf v -hi\ (2.8) 

K v = J A^-, (2.9) 

Hi 

where we have also introduced the projection operator hf along the trajectory. The 
integration over Ci($) cannot, unfortunately, be performed in an exact manner. As will 
be explained in detail in subsection 2.3, we will organize it perturbatively with respect 
to the order of t-derivatives. This expansion scheme is expected to be reliable when 
the accelerations of the D-particles are small compared with the D-particle mass scale. 
Diagramatically, it corresponds to the picture (see Fig. 2) that the dominant interaction 
occurs at a point on the worldline and the corrections due to the extendedness of the 
boundaries are subsequently taken into account by (^-integrations. 




+ £j - integration 



Fig. 2 Diagram showing our approximation scheme. 
Splitting of string coordinate 
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To perform the integration over X^(r,9), it is convenient to first split it into the 9- 
independent part Xq{t) and the rest called £ M (r, #): 

X»{r,0) = XZ(r)+e(r,9), (2.10) 
d6(*(r,0) = 0. (2.11) 

Then the action separates without surface term into 

S = ± J ±( r arX {r))* + f <Pz(d a i)* . (2.12) 

The constraints at the boundaries now read 

XS(Ri) = FM+Wtfl, (2.13) 

i"(Ri,e) = ft(ti)cm + n?(e), (2.14) 

r rl9 

where Qf >0 = / ^W) , (2.15) 

C%(0) = a?(o)-afa. (2.16) 

Note that by definition Q(9) itself has no constant part, namely, / d9Q(9) = 0. 
Boundary conditions 

We now formulate the boundary conditions that follow from the consistency under the 
variation of the action (|2.12| ). As we have to take into account the constraints ( |2.13j ) and 
( p. 14 ) already existing at the boundary, the analysis is somewhat involved. It is easy to 



see that the boundary conditions compatible with the equations of motion in the bulk 

are 

drXo(Ri) ■ 5X (Ri) = 0, (2.17) 

dri(Ri,e)-5i(Ri,e) = o. (2. is) 



As has already been stated, we will treat 0(( 2 ) terms in (|2.4j ) (which involve derivatives 
higher than the second) perturbatively. Therefore, these conditions need be imposed up 
to 0(C). As for Xq, the constraint ( |2.13|) at the boundary says that the only possible 



variation up to this order is of the form dX^Rj) = f^t^Sti, i.e. in the direction of 
the trajectory. Similarly, ( |2.14j ) dictates that (since d(9) is considered small) possible 
variation of ^(Rt, 9) is also along the trajectory. Thus, using the projection operators 
introduced in (12. 81) and (12.90, we can write the conditions as 



SXS(Bi)P^ ti = 0, d r Xg{R i )hr, ii = 0, (2.19) 
8t{R h 9)P^ = 0, d T fr(R i ,0)h mi = Q. (2.20) 
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To see the content of the conditions for X(}(Ri) in more detail, it is useful to further split 
Xq (r) into the classical part X£(r) that satisfies d 2 X^(r) = 0, X^(Ri) = ft(ti), an d the 
quantum part. The desired splitting can be written as 

X "(r) = + (2.21) 

V47T 

X&r) = - s ^mti)l^ + mh)]n r -Y (2.22) 

where s = In — . (2.23) 
a 

Putting this into (|2.19|) , we get 

= 6a%{Ri)P mi , (2.24) 
= -4-(/ 2 (*2) - h{ti)Yh mi + -i=9ra^(i2i) V,i ■ ( 2 - 25 ) 

The first of these equations simply says that x$ should satisfy the Dirichlet condition for 
the direction transverse to the trajectory at each end. The second, on the other hand, 
dictates that we should impose the Neumann condition for Xq in the tangential direction 
and at the same time demand, for consistency, 

(/ate) - fi(h))J?(h) = (/ 2 (t2) - fi(ti))J?(t 2 ) = . (2.26) 

These additional conditions may look odd at first sight, but their physical meaning is 
clear: The quantity f^ih) — fx(ti) is conjugate to the momentum transfer for the D- 
particles and the condition simply says that to 0(Q there should not be any momentum 
transfer along the trajectory so that the D-particles should be able to move freely along 
the trajectory^. We shall see later in sec. 3 that such conditions arise also from the 
requirement of BRST invariance for the boundary state representation of the interaction 
vertex. 

The remaining analysis for £ M is straightforward and we get the same type of boundary 
conditions as for Xq. So if we define the total quantum fluctuation £ M (r, 6) by 

e(r,9) = 4= + £>(r,0), (2.27) 

V47T 



we can summarize the boundary conditions as 



PrU^O) = 0, (2.28) 
K w d r URi,0) = 0. (2.29) 



J Remember that the effect of the extrinsic curvature starts at C(C 2 )- 
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Orthonormal basis 

The path integral over the quantum fluctuations £ M (r, #) will be facilitated if we set up 
a suitable orthonormal basis. In the one particle case||, there is a unique natural such 



frame, which is spanned by the unit vector u^ = f^/y—f 2 tangent to the trajectory and 
the remaining D — 1 transverse, mutually orthogonal unit vectors. Then the boundary 
conditions are neatly separated into the Neumann in the tangent and the Dirichlet in the 
transverse directions. 

For the case at hand, however, we have two ends and the Neumann directions at these 
end points are different. Thus, we may consider two possible orthonormal frames, to be 
defined below. 

Since f- 1 are expected to be time-like^, let us first introduce the time-like unit vectors 
Ui defined by 

^ = u? = -l. (2.30) 

V-hi 

In the generic case, U\ and u 2 are non-degenerate and define a plane, which we call "the 
trajectory plane". It is easy to show that the D-dimensional Lorentizan inner product 
Ui ■ u 2 takes the values in the range oo > — U\ ■ u 2 > 1, where the equality on the right 
holds if and only if the their spatial parts agree, i.e. if U\(ti) = u 2 (t 2 ). Therefore we can 
parametrize — U\ ■ u 2 in the form 

— U\ • u 2 = coshx, (2.31) 

where % i s re& l an d positive (by convention) and it depends on fi symmetrically. This 
parameter x will be of utmost importance as it will be shown, in subsection 2.4, to carry 
the essential information of the scattering. Its physical meaning will also be spelled out 
there. 

Now we can define the unit-normalized space-like vectors u 2 and Ui in the trajectory 
plane, which are orthogonal to U\ and u 2 respectively: 

u 2 = — - — (u 2 — u\ cosh y) , #2 = 1, u\ ■ u 2 = , (2.32) 
smh x 

u\ = — - — (u\ — u 2 cosh x) 1 ^i = l) u 2 -ui = Q. (2.33) 
smh x 

Either {ui,u 2 } or {u 2 ,ui} will do, but we will choose the former for our basis in the 
trajectory plane. The latter then is expressed in terms of the former in the form 

u 2 = u\ cosh x + u 2 sinh x , (2.34) 
ui = — {u\ sinh x + u 2 cosh x) ■ (2.35) 



3 When the D-particles are put on-shell after quantization, this will be automatic. 
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By appending D — 2 space-like orthonormal vectors uj which span the space transverse 
to the trajectory plane, we complete our basis as 



{e^} = { Ul ,u 2 , Ul } , A=(i,I), i = l,2, 1 = 3, 4,. ..,£>. (2.36) 

Accordingly, we will define the components of in this frame as 

U = (2.37) 

that is, £i = Mi ■ £, £ 2 = m 2 ■ £ and £j = uj ■ £. With this notation, we have 

= flA ^ = aW = aV/, (2.38) 
^ = Vab = (-+,+, •••,+), (2.39) 

for arbitrary Lorentz vectors a^,b^. 

It is now a simple matter to express the boundary conditions (|2.28 ) and (|2.29|) in our 
basis: 

d r £i(R 1 ,e) = 0, (2.40) 

URuO) = 0, (2.41) 

a(i? 2 , 0) sinh x + 6(^2, 9) coshx = , (2.42) 

0rfi(iM) coshx + d r &(R 2 ,6) sinhx, (2.43) 

URi,Q) = UR2,e) = o. (2.44) 

Note that these conditions are similar to but more involved than the ones that occur in 
the case of usual open string in a constant electromagnetic field [12, 13, H| or in the case 



of D-particles moving in the same directions ||. In the present case, the two ends couple 
to background trajectories which point in different directions and hence the boundary 
conditions for £i and £2 cannot be fully disentangled. 

2.2 Path integral over the string coordinates 
Preliminary 

We shall now perform the integration over £ M (z). Taking into account the boundary 
constraints ( |2.13| ), ( J2.14 ) , the action can be written as 



S = S cl + S d ., + S, + S al , (2.45) 
S, = ^(A -/,)!, (2.46) 
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S d - q = ^(/ 2 -/i)l(fi 2 ,o-fii,o)^ (2-47) 
a's 



S q = — ]t?z{d a (-)\ (2.48) 
S cst = -i J d9j2^(Ri,9) (e(R h e) - ft(ti)(i(6) -W(9)) , (2.49) 

i 

where (/2 — fi)± denotes the components orthogonal to the trajectory plane. In (|2.49|) 
we have introduced sources z/ M (.Rj,#), without constant parts, at the boundaries in S cs t, 
which will be integrated to produce 5- functions enforcing the constraints ( |2.14| ). (( |2.13| ) 
is already used in (|2.47|) .) 

The path integral over ^(z) then can be written in the form 

V e = J V^ e E z , (2.50) 
E e = ^ J d 2 zi(z) ■ d 2 £{z) + i J d 2 zi(z) ■ J{z) . (2.51) 

Although the source J(z), representing h>^(Ri,6), is active only on the boundaries, we 
shall, until appropriate time, take it to be a general function of z. This will allow us to 
compute the requisite Green's functions on the annulus satisfying the rather complicated 
boundary conditions imposed on £(z). Also for a while we shall suppress the component 
indices for £ and J unless needed. 

In the following, it will be convenient to switch from the variable r to a more natural 
proper-time variable p defined by 

p = m-^, 0<p<s = ln^. (2.52) 
.Til -fti 

The Laplacian and the integration measure take the form, d 2 = Ri 2 e~ 2p (d 2 + dj) and 
Jd 2 z = R 2 f 2 *d6f°dpe 2 e. 

We now expand various quantities in double Fourier series in 6 and p and then integrate 
over the modes. First expand £ and J in angular Fourier series: 

£(p,0) = £ (x n (p)+<(p))-=, (2.53) 

n V2 V27T 



i pinB 

n V2 V27T 
From the reality conditions, x n and x' n must be real and satisfy X- n = x n , x'_ n = —x' n . So 
the independent integration variables are x n and x' n for n > 1 and Xq. Then the exponent 
E% becomes 



E £ = / dp 

5 47ra' Jo H 



£ x n (-d 2 p + n 2 )x n + £ x' n (-d 2 + n 2 )x>, 



2\„/ 

n=0 n>l 
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+iRl ( S dpe 2p 
Jo 



^ X njn + "^nin 
n=0 n>l 



;2.55) 



Since the primed system is identical to the unprimed (except for n = part), we will 
hereafter exhibit the unprimed part only. 

Next we wish to make a Fourier expansion of x n (p) with respect to p in the interval 
[0, s\. Here we must take into account the boundary conditions ( 2.40| )~( [2.44 ), which in 
terms of x n (p)'s read 

3 p x llft (0) = 0, (2.56) 
rr a , ft (0)=0, (2.57) 
xi >n (s) sinh x + x 2 , n (s) cosh x = , (2.58) 
d p x 1>n (s) coshx + d p x 2 , n (s) sinhx = , (2.59) 
a*,n(0) = xi, n (s) = . (2.60) 

The expansion for xi >n , which satisfies the usual Dirichlet condition at both ends, is 
standard. But finding the appropriate Fourier expansions for Xi >n (p) and x 2i „(p), obeying 
complicated boundary conditions, is rather difficult. To circumvent this, we will use the 
following trick. The idea is to first choose some appropriate complete basis, which satisfies 
a certain boundary condition, and then realize the actual boundary conditions ( 2.56j ) ~ 
02.60T) by introducing extra source terms Aj n (p), to be discussed in detail later. For the 
moment, we include them in j n (p) and proceed. 

To find a convenient as well as consistent Fourier basis, we must take into account 
that the conditions ( |2.56|) ~ (|2.60|) are not periodic in the interval [0, s]. This suggests 
that we should extend the functions into the enlarged interval [— s, s] as even functions 
and then make the usual Fourier expansion. It reads 

x n{p) = —j= [a n ,o + ^2 a n,m cos k m p , (2.61) 

V s \ m =l 1 

TXTfl 

where k m = . (2.62) 

s 

Note that this automatically satisfies the Neumann condition at the ends of the original 
interval, i.e. at p = 0, s, which is still consistent with the true boundary conditions. We 
will use this as our base and append corrections to realize the true boundary conditions. 

So we first compute the path integral for this basic case. Substituting ( |2.61p into 
( p. 55 ), we get for the mode n (with the omission of the primed modes), 

n 2 1 00 
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+ia 2 Clnfijnfi + a n,mjn,m , (2.63) 



\ m=l 

where 



s Jo 



j n ,o = ^= I dpe 2 ?j n (p) , (2.64) 



Jn, 1 



J q dpe 2p j n (p) cos k m p. (2.65) 



The Gaussian integrations over the modes a n> o and a„ im are now trivial and it produces^ 

„ oo 

J da nfi J] da m e E ^ ~ D^(s)e E ^, (2.66) 



m=l 



where 



m-) ~ 1 fi (£ + «r" 2 = 1 ( ^i)" /2 , (2.67) 



E(j) n = -na'Rf 



1 2 °° 1 

n 2 ^ fe^ + n 2 



= -na'Rt (" dp f dp'e 2 ^ Jn (p)N n (p, p')j n (p') . (2.68) 
Jo JO 

Here, N n (p, p') is the Neumann function for the operator — d 2 + n 2 given by 

2/1 ~ cos k m p cos k m p'\ 

(coshra(s — (p + //) + coshn(s — \p — p'\)) (2.69) 



2n sinh ns 



with normalization (— + n 2 )N n (p, p') = 5(p — p'). These expressions are derived with 



p 

the aid of the formulae 



n = — -> 2 - ro 

Tn J 7ry 



m=l 

x cosmi 7r cosha(7r — 1 



= - (0<x<2vr). (2.71) 

m 2 + a 2 2a sinh an 2a 2 

We must now take the product over n. Due to the presence of the constant mode ao.o, 
the contribution for n = is divergent. As we shall see shortly, for the actual boundary 
conditions for our problem, such divergences will cancel. Thus, for now we write n = 



7 Here and hereafter, we will not be concerned with the overall constant. If desired, it can be most 
easily obtained by the open channel operator method, to be discussed in section 4. 
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as n = no, where it remains explicitly, and later check the cancellation. The front factor 
then becomes (including the contribution of x' n (p)), 



oo 

2 



n=l 



II V (i**)- 1 , 2.72 

where we used a formula n^Li sm hns = v2r](is/ir) yielding the Dedekind //-function 
rj(x) = e mx l 12 n^Li(l — e 2mnx ^ an( j employed the standard ^-function regularization, 
such as n^=i A = A~ l l 2 and n^Li n — v^vr, where appropriate. 

Transverse sector 

Let us now describe how one can realize the true boundary conditions (|2.56|) ~ (|2.60| ) by 
introducing extra source terms. 

We begin with the case of transverse components xi >n (p), which satisfy the Dirichlet 
conditions at both ends, namely xi jn (0) = x/ >n (s) = 0. In this case, we need not, of course, 
use any trick. Omitting the subscript /, we may simply expand x n (p) in the Fourier sine 
series 

f2 00 



x n{p) = \ -J2d njm smk m p, (2.73) 

» S m=l 

and go through the procedure similar to the previous case. One easily finds 

D {D \s) ~ s- x/2 -q{is/<K)- x (2.74) 



for the front determinant factor, and in the exponent of ( 2.68|) the Neumann function is 
replaced by the Dirichlet function 

D n (p,p') = — — — (- coshn(s -(p + p') + cosh n(s- \p-p'\)) . (2.75) 
In sinh ns 

This suffices for the Dirichlet case. However, in preparation for handling the more difficult 
case of £1-^2 system, it is important to check if the same result can be reproduced by our 
more versatile trick, to be described below. 

In terms of the modes defined in ( |2.61| ), the Dirichlet conditions can be implemented 
by the insertion of the following integrals representing £-f unctions 8 : 

J dp n exp f ^=(a„ i0 + V2 ]T a n>m ) j , (2.76) 



m=l 

8 One might worry that imposition of the Dirichlet conditions on top of the Neumann conditions 
already built in is overconstraining. Classically, it is a legitimate concern. But in the path integral 
calculation, the effect of the extra Neumann conditions is to restrict the behavior of the functions only in 
the infinitesimal vicinity of the ends of the interval and this modification is of measure zero in the space 
of functions to be path-integrated. 
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dp' n exp f^Ko + y/2 £ a n , m (-l) m )"j . (2.77) 
\v s m =i / 

The terms in the exponent can then be regarded as an addition of extra source terms of 
the form 

Aj n ,0 = -^(/i n + /4J, (2.78) 

a 

Aj„, m = -^=( M n+(-l)VJ, (2-79) 
which, in the light of definitions (|2.64|) , (|2.65|) , can be written simply as 

Ri 

In other words, we are putting random sources at the ends to enforce the Dirichlet con- 
ditions. Replacing j n (p) in ( |2.68| ) by j n (p) + Aj n (p), we get, after some calculations, the 
following additional terms for the exponent: 

AE(j) n = =-^a'{Nyi + ^—n' 2 \ 



= —^(pJ(p)+p'J(s-p)). (2.80) 



TV 



<. N o> 



N° ( f i x 2\ 

+wa'Ri , J (J n (s)) 2 tanh 2 ns + (J n (0) — J n (s) 

sinri tis \ \ cosh t\s 



(2.81) 



where 



N° ee N n (0,0) = N n (s,s) = ^^, (2.82) 

n sinn ns 

J n {s) = f S dpe 2p Jn {p) cosh n{s-p), (2.83) 
Jo 

J n (0) = I* dpe 2(, j n (p) cosh np. (2.84) 
Jo 

Modified sources p n and p! n are introduced in perfecting the squares but the Jacobian fac- 
tor is unity for this rewriting. Upon integration over these sources, the determinant factors 
cancel, except for a factor of n. Taking the product over n, we then get n (n^Li n ) 2 ~ n o, 
up to a constant. This extra factor of n cancels l/n in fl2.72p and we reproduce fl2.74p . 
On the other hand, the remaining terms in ( |2.81| ) give, after some rearrangements, 

- ira'Rt (- f dp f dp'e 2 ^j n (p) Jn (p') 2 costing - (p + p'))) , (2.85) 
V Jo Jo 2n sinn ns / 

which when added to the original piece fl2.68p precisely effects the change N n (p, p') — > 
D n (p,p'). This demonstrates that our extra-source method indeed works nicely. 
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£i-&2 sector 

With this warm up, let us now proceed to the main task of incorporating the boundary 
conditions for the £1-^2 system. Since the basic method has already been described in 
detail, we shall only give the essence. 

The extra source terms for xi >n (p) and x 2j „(p) to be added now take the form 



Aii(p) = ^Ms-p)smh X , (2.86) 



^fn(p) = %^(fi5(p) +fi' n 5(s-p) coshx). (2.87) 

Substituting the shifted expressions j l n (p) + Aj^(p) into the exponent ( |2.68|) , where now 
we must replace j n (p)N(p, p')j n (p') by jn(p)r)ij N (P, P')jl(p') witn tlie metric ^ = (-, +), 
we get after some calculations 



E n (j) = -Aira'(N° n j2l + ^§^p' 
\ cosh ns 

+ixa'R\ 



■ r 1 r J:W*(s) 
n sinn ns cosh ns 

cosh ns ( ~ , , , , . . , cosh v . . 

+ — — — — Jl cosh X - j; 1 sinh X r^-Jnis) 

nH„ smh ns \ cosh ns 



2- 

,(2- 



where the quantity if n is defined by 



H n = sinh 2 ns - sinh 2 X = \e 2ns {l - e 2x " 2ns )(l - e ~ 2x - 2ns ) , (2.89) 

and Jn(s), J7^(0) are defined just as in ( |2.83| ) ~ ( |2.84j ) with f n (p) in place of j n (p). Inte- 
grating over the sources p n ,p' n , and forming the appropriate products over n (with the 
contribution of the primed system as well), the determinant factor becomes 

£>(i2)( s ) ^rjfis/rfe 2 */ 12 TT (1 — e 2 *- 2 -)- 1 ^ - e -2x-a»-)-i , ( 2 .90) 

smh x n=i 

where the factor n^s/ sinhx is due to the n = mode. This must be multiplied by the 
already existing factor D^ N \s) 2 ~ {l/n 2 ) s)n{is/ , K)~ 2 for two degrees of freedom. The 
factor n^s and the ^-functions cancel and we get 

1 00 

£,(12) ( s ) „ _^_ e 2 S /12 TT (l _ e 2 X -2n S) -l (l _ e -2 X -2n S) -l ,3 gl) 

smh x 

As for the Green's function, we can read off from Q2.88 ) the corrections AGj 3jTl to be added 
to i]ijN n (p, p'). Omitting the details, the final result for the full Green's functions is 
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Gn. n (p,p') = — (sinhn(p + p') + sinhn|p - p'\) 
In 



sinh ns cosh ns . . . ,. . . ... . . 

— (cosh n(p + p ) + cosh n(p — p )) , (2.92) 



2nH n 

G\2,n{p,p') = (sinh n(p + p') - sinh n(p - p')) , (2.93) 

G 21>n (p,p') = S -^^(smhn(p + p') + smhn(p-p')), (2.94) 

G 2 2,n{p,p') = g - (sinh n(p + p') -sinh n|p-p'|) 

sinh ns cosh ns . ,. . ... , , 

(cosh n (p + p ) — cosh n(p — p )) . (2.95) 



2n#„ 

One can check that they satisfy (— <9 2 + n 2 )Gij tTl (p, p') = %-5(p — p') and that in the 
limit x they reduce to the expected forms, namely, Gii >n — > —N n , Gi2,n —* D n , and 
G\2, n -, Gi\, n ~~ * 0. Moreover, they are of such forms that we can take the limit n — > to 
obtain the Green's functions for the zero mode as well. 

Summary of the lowest order amplitude 

Let us now assemble the results so far obtained and write down the 0-th-order amplitude, 
i.e. yet without the corrections at the boundaries. It consists of the contribution from 
the classical part e~ Scl with S c i given in (|2.46| ), the one from each transverse coordinate 
£/ given in ( |2.74j ), and the one due to the £1-^2 system, ( |2.91| ). To these we must add 



the contribution of the b-c ghosts, which is , as usual, D^ 9h \s) = n{is/n) 2 . Finally, we 
must integrate over the proper time s. Since we are dealing with a cylinder, the proper 
measure should simply be J °° ds. Thus, using the explicit form of 77 (is/71"), we get 



Vo(A,/ 2 ) ~ / dse 2s e( D ~ 26 >/ 12 s-( D - 2 V 2 e-V 2 -VV 2a ' s 
Jo 

•-L- fl (1 - e- 2 -y iD ^ (l - e-^)" 1 (l - e-— ) ^2.96) 



71=1 



as our lowest order approximation to the amplitude with fixed trajectories. In terms of 
the usual 6- and ^-functions, it can be written compactly as 

V (/i,/ 2 ) ~ rds S -^- 2 )/ 2 e -^^)i/ 2 ^(r)-^- 2 )^j4, (2.97) 

JO vi[v\T) 

where r = is/rc and v = —ixj^- 
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2.3 Corrections at the boundaries 



We shall now compute the corrections due to the extendedness of the boundaries, i.e. to 
the fluctuations of the geodesic coordinates 0(0)- 

The first task is to derive the effective boundary action for d(9) arising from the 
integration over the special sources i/ M (.Rj,0) that we introduced in (|2.49|) to enforce the 
constraints ( |2.14| ). These sources have until now been represented by the general source 
J A (z) (or its Fourier components j£) and the integration over the string coordinates 
^(z) produced an exponential factor with the exponent quadratic in J^'s connected by 
the Green's functions. 

To apply this to the present case, we must recall the following: (i) ^(Ri, 0)'s do 
not have constant parts, since they are designed to couple only to the non-constant part 
of (ii) Because, at each boundary, components of £ M orthogonal to the direction of 
the trajectory vanish by the boundary conditions, the only components of v^{Ri,9) that 
actually couple to £ M are the ones along the trajectory, namely u(Ri, 9)-U\ and u(R2, 9)-U2- 
The point (i) says that we only need a Green's function in the space of non-zero modes 
(with respect to 9). Furthermore, the point (ii) tells us that only the Green's functions 
with components in the trajectory plane will be required. 

Thus we form the following two-dimensional Green's function on the worldsheet with- 
out the zero mode: 

G l3 (z,z') = G ij (p,9;p',9') = £G 0> (/>,p) 008 71(0 -00, (2.98) 

n>l 

where GV, >n (p, p'Ys are as given in (|2.92|) ~ Q2.95Q . This satisfies —d 2 Gij(z, z') = n5 2 (z, z'), 
where 5 2 (z,z') = (1/ 'R\)e- 2p 5{p- p'){5(9 - 9') - (1/2tt)) is the delta function in the space 
of non-zero modes. Then defining 

z/(0) ee u(Ri,9)-u\ (2.99) 

one finds that the exponential factor quadratic in the source takes the form 

exp J dpdp'd9d9'j i (p, 9)G ij (p, 9; p', 9')j j (p', 6')^ , (2.100) 

where 

j\p,9) = v\9)5(p)+v 2 (9)6(s- p)cosh X , (2.101) 

f(p,9) = v 2 (9)5(s-p)smh X . (2.102) 

Using the explicit forms of Gy >n (p, p'), this exponential can be written as 

exp(-a' J d9d9'u l {9)g iJ {9,9')u J {9'^j , (2.103) 
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with 

r, m nt\ ^ sinhns / coshns coshv \ „ A ln 1ri ,. 

Note that at this stage the Green's function, Gij(0, 0'), has become completely symmetric 
in This must be the case because the indices % — 1, 2 for z/(#) now refer to the two 

trajectories symmetrically. 

We must further add to the exponent the remaining terms linear in the sources that 
directly couple to G (#)' s ( see (|2.49| )). It is convenient to separate out the parts containing 
v % {&) and write them in the form 

ifde\J2 - n(9)j , (2.105) 

where 

z t (0) = fT^m + ^-0=cm 3 + ■ ■ ■ , (2.io6) 

Sl(0) = Y.v^ d ) p r(^K hV Cm 2 + --) ■ (2-107) 

(2.108) 

In the above, d(0) n means that the constant part contained in d(6) n is removed. Then 
integration over u l (6) yields 

(det^)- 1 / 2 exp (^-j J dBdB'Z i {B)Q- 1 {d,d') ii Z i {d') \ , (2.109) 

where the inverse Green's function Q^ 1 {9^9 1 ) is easily obtained from (|2.104|) as 



g-\e,e>) = -If^^f cosh ^ ~ co ^ ) C0S n(9 -e>) . (2.110) 

7r 2 ^ nH n \ — cosh x cosh ns J 

We have now completed the calculation of the effective interaction of Q(j9)'s. It consists 
of the factors e~ s ^ from ( ggTp , e^/^W from (|2~T05l) , and the factor ( ^T09| ) just 



computed. Note that so far no approximation has been made. 

To proceed further, we now need to make an approximation. As was already an- 
nounced, we will consider the case where the accelerations fa ( and higher tj-derivatives 
of fi{ti) ) are small and treat the terms containing them as perturbations. Specifically, 
in this article we will keep only the terms linear in the extrinsic curvature = f v . 
Then, the expression we need to deal with is 
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V corr = (det G)- 1 / 2 exp J dOdO'^i (,(6)0-^6,6^^(3(6') 

1 f d6 



+^E/ d6u fl (R l ,6)PrK i ^(6) 2 



(2.111) 



The Gaussian integration over d(6) yields the following effects. First, it produces the 
determinant factor, which consists of (det(7) 1//2 which cancels (det^) -1 / 2 in front and 



(n(-v«r l/2 ) 2 (n(-v«r l/2 ) 2 = f^fv^, (2.112) 

n=l n=l 

which is obtained by the Fourier mode expansion and the ^-function regularization. The 



factor \J—hi will be important in changing dt, into the proper time dr, ~ dt^—hi. Sec- 
ond, the terms quadratic in d(6) are replaced by their quantum averages. The relevant 
correlation function {0(6)(j(6')) is given by 

In 1 

<C<(W)> = - = r— GaW). (2.113) 

In particular, using the explicit form of the Green's function, the ones at the coincident 
point are computed as 

r\ / q / / 00 1 \ 

UO) 2 ) = — Gu(6,6) = ^- lne-£ — -(cosh2 X - e" 2 -) , (2.114) 

where we have regularized the divergent sum J2 n >i(^-/ n ) ~^ En>i( e_£ 7 n ) = — me - Since 
this is independent of 6, we immediately have (0(d) 2 ) = 0, so that the last term in (|2.111|) 
does not contribute to this order. As for the term containing f2 — fi, we get a divergent 
term as well as a fininte correction. The former is of the form 

-i(/ 2 -/,k-f^-^Vn<. (2.H5) 



s \ hi h\ 

If we remember that we have, in the lowest order amplitude, the classical part of the form 
exp(— (l/2a's)(/2 — fi)±), we see that ( |2.115| ) can be absorbed by the renormalization of 
the trajectories (subscript R stands for renormalized quantities) 



ft = + (2.H6) 



$ft,R = -a'-^lne. (2.117) 
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It is gratifying that ( 2.117 ) is exactly the form found for the one D-particle case treated 
in |B 

Putting all the findings together, we now have the final form of the amplitude for fixed 
trajectories, valid up to the first order in fi\ 

V(/i,/ 2 ) = /n^V^V^/f dse 23 ^ - 2 ^' 1 ^-^- 2 ' 2 ^-^-^ 2 *'* 



[] l-e" 2ns ' ' (l-e- 2ns+2 x) 1-e 



sinhx^Ji 

. e (/ 2 -A)x-AK ? ^ 2118 ^ 
1 / \ 00 1 

where AK* = - - -±- V — -(cosh 2 X - e' 2ns ) . (2.119) 
2.4 Amplitude for quantized D-particles 
Quantization of the collective coordinates 

Having obtained the amplitude for fixed trajectories, we now quantize the trajectories 
themselves. The basic formalism developed in || works with a slight modification and 
hence we shall only describe the essence and refer the reader to for further techinical 
details. 

To the accuracy of our approximation, the action for a relativistic D-particle is of 



the form So = —im Jq 1 dty — f 2 ) (with m the bare mass), which, following [|15|| , can be 
effectively turned into a more tractable quadratic action 

S = f drhf 2 {r)-m 2 ). (2.120) 
Jo I 



Here, m is the renormalized mass and r(t) = (1/m) / * dt' J —h(t') is the (rescaled) proper 



time. Then the t-derivatives can be rewritten into r-derivatives in the manner f{t) / J—h(t) 
= f(r)/m, f(t)/(—h(t)) = f{r)/m 2 , etc. Hereafter, dots refer to r-derivatives. The 
quantum amplitude is then obtained by A = Jq 00 dT J T>f{r) exp(iS)V (/), where V (/) is 
the amplitude ( or the "vertex operator" rather) for fixed trajectory. 

For the present case, V (/) computed in the previous subsection depends on /j and fa 
as well as on f { . To separate out the dependence on the derivatives, it is convenient to 
replace them by new variables i>„i> where n{= 1,2) refers to the number of derivatives, 



We have also confirmed that divergences that occur at the second order in are precisely absorbed 



by the renormalization of \f—h\ \f—hi with the same prescription (2.117). This is a highly non-trivial 
check of the consistency of our treatment. 
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via the insertion of unity: 



1 = J]J dv n,i J]Jdu nA expiiJ2^n,i( v n,i- dn fi( T i))) ■ (2.121) 

n,i n,i \ n,i J 

Then the path integral over fi{j{) we need to consider at the first stage is 

^ = /n^)ex P fiE^ T< ^(/^)-^ 2 )) 

• J dndT 2 exv (-(/ 2 (r 2 ) - h{r l ))H{2a's) + (/ 2 - f,) ± ■ AK(v n ,)) 

•exp \-iJ2 UJ n,i-d n f i {r i ))) . (2.122) 



Let us denote by fa and // the initial and the final positions of the D-particles respectively, 
and decompose /f (r,) into the classical part /^(Tj) and the quantum part ft( T i)'- 



where 2/i 

7,(o) 



n + ft, 



Hi 

Ti 

fi fit 

fi(Ti) = 0. 



(2.123) 
(2.124) 
(2.125) 
(2.126) 



Then due to the boundary condition ( |2.126|) the classical and the quantum parts separate 
in the action: 



- / dr(ff - m 2 ) = A. - -m 2 T t + - drf, . 
2 Jo KJl ' 2Ti 2 2 Jo 1 



(2.127) 



On the other hand, for the term (^(t^) — /i(tl))j_ such a separation appears difficult at 
first sight. This problem,however, is solved by going to the momentum representation. 

Let pi and p\ be the incoming and the outgoing momenta of the D-particles respectively 
as in Fig. 3. 



P 2 



Pi 



P 2 



Fig. 3 Scattering diagram in momentum space. 
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Then the factor to be multiplied for going to the momentum representation is 

expOi ■ fi - ip[ ■ f[) exp(ip 2 ■ f 2 - ip' 2 • f 2 ) 

= exp(-i( yi ■ pi + y 2 ■ p' 2 )) exp(i(/ • {pi - p[) + f 2 ■ (pi + p 2 - p[ - p 2 ))) , 

(2.128) 

where we have defined / = j\ — f 2 . Integration over f 2 is now trivial and produces the 
momentum conserving (^-function 5(pi +p 2 — p[ —p' 2 )- As for the integration over /, it is 
of the form 

J = f d D f e if-(pi-p'i) e -(f+9)l/(2a's)+(f+g) x -AK ^ (2.129) 



where g = - f^r 2 + fi(n) - f 2 (r 2 ) . (2.130) 

So by shifting / + g — > /, we get 

= e -ifl-(pi-Pl) /" d D^ e ifi Pl - P ' 1 ) e -fl/(2a's)+f ± -AK 

~ s^-^e-^^-pDe-^'^CPi-Pi-^)!^^^ .p^,). ( 2 .131) 

^ 2 ((K — Pi)||) is the two-dimensional (^-function for the components in the plane spanned 
by the vectors f^i's. Note that since g in the exponent is linear in y^, the classical and 
the quantum parts are now separated. Note also that the front factor cancels s _ ( Z)_2 )/ 2 
already present in the amplitude. 

Next we perform the integration over y^. Gathering together the terms containing y^, 
the integral is 

I y = J d Vl dy 2 exp jz ^ ^— y i yi Jj 

~ {T l T 2 ) D ' 2 ex V i^-Y,^^+p[T l + T l {p l -p[)f^ , (2.132) 

where we have used momentum conservation. 

Next we deal with the quantum fluctuations fi. The integral to be performed is 

rT -J2 



J £*LJfi(0)=fi(T)=0 \2Jo V / 

(2.133) 



This is precisely of the type that occured in || and can easily be evaluated by using the 
one-dimensional Green's function satisfying the Dirichlet conditions at the ends of the 
interval [0, T]. As it was shown there, with the standard (^-function regularization, the 
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contributions of d n f{ji) average out to zero for n > 2. Skipping the details, the result is 
(using again the momentum conservation) 

t rmrrs-Dn .(Pi-Pi) f Ta i 2 \ 

I - f = {T X T 2 ) 1 exp I 2^ 1 gjT ~ z r \ l ~~ ~2J 2X ^ Pi ~ Pi > ^ ~ " j ' 

(2.134) 

Combining ( |2.132| ) and ( 2.134j) , the front factors cancel and the exponent becomes 



E (-~"i,i(Pi+P'd + \(P' 2 1-Pl)r* - \v'-T^ . (2.135) 

If we now put back the factors exp(z ^2 n iUJ n ,i v n,i) ( see Q2.121D ) and integrate over u n ^, we 
get the following 5-functions: 

This means that effectively we can replace fi by \(j>i + p^) and set to zero, exactly 
the same rule that we established in [Bl. Therefore the correction AK that remained in 



( |2. 13 1| ) actually vanishes. Another consequence is that what has been referred to as the 



"trajectory plane" should now be understood as the one spanned by the mean momenta 
\{p'i +Pi) and \{p' 2 +P2). Then when the D-particles are put on shell, we easily see that 
the momentum transfer k = p[ — pi = — (jp 2 — P2) is orthogonal to this plane since 

k-^(p' 1+Pl ) = I(^ + m 2 )- I( p 2 + m 2 ) = 0) (213?) 

-k-\(p' 2 + P2) = l -{p'l + m 2 )- l -{pl + m 2 ) = V. (2.138) 

This is quite consistent with the presence of S 2 ((p[ — Pi)\\) in ( |2.131| ) and allows us to 
replace (p[ — Pi)\ in the exponent by (p[ — p\) 2 ■ 

The remaining integrations over the interaction points Tj and over Tj are trivial and 
they produce the four propagator legsQ, to be removed from the proper amplitude. 

We may now write down the quantized proper amplitude in complete form. Omitting 
the momentum conserving ^-function and with the understanding that the D-particles are 
put on shell, it reads 



AUpupM*) ~ / dse 2s e^- 2 ^l l2 e~^ s l 2 ^-^ 2 

Jo 

ft (1 - e- 2 -r (D - 4) (1 - e-^r 1 (l 

sinhx - ' v ' v 



n=l 



(2.139) 
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where \ is defined by 



coshx = —u\-U2, (2.140) 
ik = y R + PS ■ (2-141) 



This, however, is not quite the correct answer for fully quantum amplitude: We have not 
yet taken into account the quantum indistinguishability of D-particles. It is clearly seen 
when we express cosh \ defined above in terms of the velocity v and the scattering angle 
8cm in the center of mass frame of the D-particles. It takes the form 

cosh , = \ + ^X CM ?l < 2 - i42 > 

1 - v 2 cos 2 {9cm/2) 

One sees that , for any v, x vanishes for 9cm = 7r, i-e. for back-scattering but it does 
not vanish for the forward scattering. For indistinguishable particles this is obviously 
incorrect. The cure of course is to add the amplitude in which Qqm is replaced by 7t — 9cm- 
This replacement must also be done for the exponential factor exp (—(a's/2)(p' 1 — Pi) 2 ), 
which in the center of mass frame takes the form exp (— 2a' p 2 sin 2 (9cm , with p the 
maginitude of the spatial momentum. More generally the replacement to be made is 
Pi P21 so that the desired quantum amplitude is 

AbosijPhP'i) = A bos (p 1 ,p2;p 1 ,p 2 ) + A bos (p 1 ,p 2 ;P2,p 1 ) • (2.143) 
Forward scattering in the infinite mass limit 



Let us check our result against the known expression ||, [16] in the special case of for- 
ward scattering for infinitely heavy D-particles. In this case, the particles are treated as 
distinguishable backgrounds and hence we should use ( |2.139p ~ ( |2.141| ). 



To make the comparison, we first turn the amplitude into the impact parameter repre- 
sentation by introducing a transverse vector and performing a Fourier transformation 
with respect to the momentum transfer k± = (p[ —pi)±. Since x does not depend on the 
transverse components, this transformation simply produces the factor 

J d D - 2 k ± e^ b e-^' S / 2 ^ ~ s (D-2)/2 e - b y { 2a>s) Q ^ 

Let us now assume that the D-particles are moving in a common direction, say 1, with 
velocities V\ and V 2 - Then the 4-momenta of the particles are given by = (£";, EjYi, 0), 
where E t is the energy of the particle #z, and they do not change as they scatter. Hence, 
the unit-normalized vectors Ui take a simple form 

Ui = IJM. (2.145) 

" - V 2 
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From this sinh \ is easily computed to be 



\Vi-Vi 



sinh x = V(m-«2) 2 -i= r-— r- — =, (2-146) 



which is exactly the expression that appears in || [16| . Applying ( 2.144 ) and Q2.146j ) to 
( j2.139| ), we recover the known result in this special case. 

3 Boundary State Approach 

In this section, we shall demonstrate that the lowest order amplitude for fixed trajectories 
obtained in ( |2.96|) can be reproduced by using the boundary state representation developed 



in [T(| To avoid possible confusion, we wish to emphasize at the outset that the 
"boundary state" that appears below does not correspond to a definite state of a D- 
particle. Rather it should properly be understood as a representation of the DO-DO-string 
interaction vertex describing the transition of a state of a D-particle into another state^]. 
This feature explicitly appears as the presence of the derivative fi, which, when the D- 
particles are quantized as in the previous section, turns into \(j>i + p'i) containing both 
the initial and the final momentum. 

3.1 Boundary state representation of the vertex 

The boundary ket \B; f) associated with a trajectory f^ l (t) is of the structure 

\B;f) = \B z ;f)®\B nz ;f)®\B gh ), (3.1) 

where the subscripts "z" , "nz" and "gh" stand for zero mode, non-zero modes and ghosts, 
respectively. The explicit forms for the non-zero modes and the ghost parts for our 
problem are given by 

\B nz] h) = exp(-^-a^(A)a_ fl ,J |0> , (3.2) 
(B nz ;f 2 \ = (0\exp(-^±aZD;(f 2 )a nt )j , (3.3) 

oo 



\B gh ) = exp \ J2 [ c -nb- n + c^ n b^ n J (c + C )| ||) , (3.4) 

(B gh \ = (TT I (fro - b~o) exp [c n b n + c~b n , (3.5) 
where D^f,) = (h^ - {Pfi . (3.6) 



10 This feature is obscured in the case of infinitely heavy D-particles, since then their states do not 
change by the interaction and the "vertex" looks like an ordinary state. 
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The matrix D v (fj) has eigenvalue +1 along /, (Neumann) direction and —1 for the trans- 
verse (Dirichlet) directions. 

As for the zero-mode part \B z ;f), it should essentially be the position eigenstate 
\f) since the trajectory is given fixed. However, as was discussed in section 2.1, we 
must also incorporate the requirement that the momentum transfer in the tangential 
direction should vanish. This turned out to follow automatically from the requirement 
of BRST invariance. This invariance demands [|T2]| L n \B z ; /) ® \B nz ; f) = L_ n \B z ; f) <g> 
\B nz ; f), where L n and L n are the closed string Virasoro generators in the usual notation. 
Then, since a.£\B nz ; f) = —D^(f)a"_ n \B nz ; f) holds, one easily finds that we must have 
p^D v (f) = —p u , or equivalently, p^u^ = on \B Z ; f), with the unit vector along f 1 . 
Thus the correct form of the zero- mode part should be, up to a constant, 

\B z -h) = Sip-uJlh), (3.7) 
(B z ;f 2 \ = (fM P -u 2 ). (3.8) 

One can check that in the case of infinitely heavy D-particles moving parallel to each 



other the ket (B.7I) reduces to the one constructed in 16 



3.2 Calculation of the amplitude 

The amplitude is now given by 



V (/i,/ 2 ) = (B-J 2 \—-\ -U Q \B;h) 



Lo + Lq — 2 

/ ds{B-f 2 \e- s ^ +LQ -^U Q \B-f l ), (3.9) 
Jo 



where 



a' 



L + L - 2 = — p 2 + ( a -n ■ OC n + a- n ■ 6i n ) 
1 n=l 

—nC-n ~r 0— n C n ) . 2 (3.10) 

n 

is the usual closed string Hamiltonian and Uq is a ghost zero mode insertion, to be 
discussed below. 

Let us now compute the contribution to the amplitude from each sector. First, for the 
zero mode sector, we have, ignoring the overall constant, 

V z = (f2\Kvu2)^- a ' s/2)p2 5{p-u l )\f l ) 

' d D pe ip ^ h)+{ ' a ' s/2)p2 5(p ■ u x )5{p ■ u 2 ) 

s -(D-2)/2 



sinh x 



e -(/ 2 -A)i/(2a' fl ) (3n) 
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The exponent only contains (/2 — fi)±, the components transverse to the trajectory plane, 
and this agrees with the discussion in section 2.1. The factor 1/ sinhx comes from the 
momentum integral in the trajectory plane, J d 2 p5(p ■ u\)8(p ■ 112), which gives the inverse 
of the area spanned by the non-orthogonal unit vectors u\ and The expression (|3.11|) 
agrees completely with the corresponding zero mode contribution computed by the path 
integral method. 

Next, we turn to the non-zero modes. Propagation through the proper time s gives a 
factor e~ 2ns inside the sum in the exponent of the boundary ket and hence what we need 
to compute is 

" ex P E le-^a^D^a^ |0> . (3.12) 

This can be evaluated by normal ordering of the oscillators. In the Appendix, we de- 
rive a complete normal-ordering formula^] by extending the method developed in fl7|| . 
Let (ai,a\) and (<2j,aj) be two independet sets of oscillators satisfying the usual (anti- 



) commutation relations 



5ij, A and B be arbitrary matrices, and 



write aAa = J2i 4 caAijCij etc.. Then we have 



exp (aAa) exp [a) Ba)^j 

= [det(l - eBA)r e exp (a\l - eBA)- 1 Ba ] ) exp(aA(l - eBA)a) 

■ exp (-a f ln(l - eB T A T )a - d) ln(l - eBA)a) , (3.13) 



where e = +1 (—1) for commuting (anti-commuting) oscillators. Setting 



C^n,fii a n,u ■ 1 — a -n,At i (3-14) 



i * 1 

n 

A ni i,iv = Dy,v{h)$nf., B n ^e u = e 2ns D fiu (fi)5 n £ , (3.15) 

we readily obtain 

V nz = II [det (l - D.Die- 2 ™)} , (3.16) 



n=l 



11 For the present purpose, we only need the formula for the vacuum expectation value, which is certainly 
well-known. The full formula given below, however, is needed in more involved computations (such as 
those of Green's functions, etc.). Since, to our knowledge, it has not been recorded in the literature, we 
provide a derivation in the Appendix. 
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where the determinant here is over the space-time indices only and we have written D { 
for D u (fj). It is easy to show that in our orthonormal basis {u\, u 2 , «/} defined in fl2.36|) . 
the Di's have the following representation: 

A = di®(-l D - 2 ), (3.17) 
dx = ( J \ y (3.18) 

d = ( cosh2x sinh2x \ , g ^ 

I — sinh 2x — cosh 1\ J ' 

where 1_d- 2 stands for the unit matrix in the D — 2 dimensional space spanned by Ui. 
Therefore 

n r t i cos1i2y sinh2v \ /r . 

= <M 2 © = sinh2 ^ cosh2 ^ J © l - 2 • (3.20) 

Note that the matrix d\d 2 is of the form of a Lorentz boost in two dimensions and hence 
its eigenvalues are e ±2x . Therefore we get 



v nz = n [ det (i - ^i^ 2 e- 2 ' 

ra=l 

oo 

(l - e - 2ns+2x ) (l - e - 2ns - 2 *) (1 - e -2-)-P-2) . (3.21) 



n=l 



Finally, consider the ghost contribution. As is well-known, due to the existance on 
the cylinder of one Teichmuller parameter and one conformal symmetry (rotation), we 
need to insert a zero mode Uq = (cq — co)(&o + b) in order to get a non- vanishing inner 
product |12|. Then, applying ( |3.13| ) with e = — 1, we get the familiar answer 



V gh = (B gh \e-^o +L ^U \B gh ) = l[(l-e- 2 n 2 - (3-22) 

n=l 

Assembling altogether we thus obtain 

V (/l,/ 2 ) ~ /' 00 rfse 2s s -(D- 2 )/ 2e -(/ 1 ( tl )-/ 2 fe))i/2a' S 
J 

g (1 _ e -2nsy(D-4) ^ _ e - 2 n, +2x ^ ^ _ e - 2 ns- 2x ) ^ 33) 



sinhx^Jl 



which, for D = 26, is identical to the lowest order result ( |2.96P obtained by the path 
integral method. 



We can also compute the Green's function Gij(z, z') (see ( |2.98| )), satisfying the rather 



complicated boundary conditions, by using the boundary state representation. Since the 
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one we need is in the space of non-zero modes, the relevant quantity is the non-zero mode 
part of the amplitude with two closed string tachyon vertices inserted 

A nz (k,k>) = (S n ,;/ 2 |e*^ (+) We ife ^ C - ) W e ife '^ (+) ^) e ifc ^ ( - ) ^| jBriz;/l)) ( 3 .24) 

where X < - ± - ) are the non-zero mode parts of the closed string cooridnate given by X^'(z) = 
-R Y^=i(l/n)(a^ n z + a Tn z ±n ). To compute this amplitude, the normal-ordering for- 
mula ( |3.13D is indispensable. As we wish to check only the non-trivial components 
Gij(z,z f ) with i,j in the trajectory plane, we restrict the momenta k and k' to have 
components only in this plane. Then, combining with appropriate use of coherent state 
method, we obtain, after a long calculation, 

A nz {k,k') = Dexp^G^z,^)^) , (3.25) 



where D is the determinant factor G3.16Q , k % stands for k-u % etc., and Gij(z, z') is precisely 



the one that we computed in the previous section by the path integral method. 

4 Operator Formalism in Open String Channel 
4.1 Amplitude in the open string channel 

In the previous two sections, we have computed the DO-DO scattering amplitude essentially 
from the closed string channel. An important question is to see how, in the low energy 
domain, the result can be understood in terms of a spontaneously broken effective gauge 
theory Q, which is based on the open channel picture. Although there is no doubt that 
this gauge theory description is basically correct and is extremely useful, how it should 
work for D-particles with finite mass (and hence with recoil) is still a non-trivial issue. As 
was already emphasized in the introduction, it must be intimately related with the "non- 
commutative nature of spacetime" and clarifying this connection would help understand 
the physical meaning of this intriguing concept. This makes the following re-derivation 
of the amplitude in the open channel a significant exercise. 

Before taking up this task, let us first perform the modular transformation of the 
lowest order amplitude ( |2 . 9 7| ) to see what to expect. Introducing the usual open channel 
modular parameter w as 

w = e 2 ^, Z^ = Z^, (4.1) 

7r In w 

and using the standard formulas for the modular transformation of the ^-functions, the 
expression ( |2.97| ) for D = 26 is transformed into 
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1 dw 

fH 



V (/l,/ 2 ) ~ ^(/x-/ 2 )i/(^a') w ,(l-)/ 2/(w) -24 



W 2 



n~ =1 (i-^ m ^)(i-w m ^) 



(4.2) 



where 



■ X 

v = —%— 



f(w) = n (1 - ™ m ) • (4-3) 

m=l 

This indicates that apparently the open string spectrum is a peculiar one in that (i) 
the Casimir energy v{l — v)/2 is complex, (ii) two types of non-zero modes exist with 
complex energy levels m + v and m — u, and (iii) an additional excitation appears with 
pure imaginary energy v. By performing the quantization in the open string channel, we 
should be able to understand the meaning of these modes and reproduce the amplitude 
in the form of the partition function 

Z open = -~Trln(L -l) = i — ffL^Trw 1 ^ 1 . (4.4) 
2 2 w \mw/ 

4.2 Quantization, spectrum and the amplitude 

Let us now perform the quantization in the open string channel 0. We will employ the 
usual strip coordinate (a, r), where < a < n. As in the calculation in the closed channel, 
we split the string coordinate into the classical part X£ and the quantum fluctuation part 
In the present coordinate system, the classical solution takes the form X^(cr, r) = 
^(f^ih) — fi(ti)) + fi(ti), which describes a string stretched from f\(ti) to f^ih)- 
Remembering that the components of j'2 — f\ in the trajectory plane must vanish by 
consistency, the energy associated with this stretching is H d = (/i(ti) — /2(^2))i/(4vr 2 a'), 
which immediately gives the factor y/.ft~/ 2 )i/( 47r2Q ') appearing in (|4.2j ). 

As for the fluctuation part, we will continue to use the orthonormal frame defined 
in (|2.36|) . Since the D — 2 components transverse to the trajectories satisfy the usual 



Dirichlet conditions and their quantization is standard, we will concentrate on the tangent 
components £1 and £2- By now the appropriate boundary conditions shoud be familiar: 

cUi(0,r) = 0, £ 2 (0,r) = 0, (4.5) 
£i(7r, r)sinhx + 62(7^) cosh x = q, (4.6) 

(7t,t) cosh x + <9o-£ 2 (tt,t) sinh% = 0. (4.7) 

(4.8) 



12 The calculation will turn out to be quite similar to the one perfomed in g (see also [lj]) but we 
will discuss the subtleties involved and their physical interpretation in much more detail. 
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It is not difficult to find the solutions of the free wave equation satisfying these conditions 
and expand £1 and £2 i n terms of these modes. A convenient way of writing the expansions 
is 



6 



UxMrp + e ~(xMrp\ CQsh * a 



E h{< 



n>l,e=± n 



(i tX " T - al n e 4A ™ r ) cos X £ n a 



A, 







+iv a 



Ux/i&rp + e -<X/*)rp\ sinh 

7 E 



X 



-a 



7X 

a e ^e* A ™ T ) sin Ma 



(4.9) 

(4.10) 
(4.11) 



where A^ are given by 



n>l,e=± n 



A; 



, -X 
m±i— . 

IX 



(4.12) 



As was already anticipated from the modular transformed form (|4.2|) , the eigenfrequen- 
cies are complex and the modes grow or diminish as functions of r. Technically, this 
comes about because the boundary conditions involve hyperbolic functions of x while 
the oscillatory part is trigonometric. Physical origin will be discussed after we finish the 
quantization. 

Since the creation and the annihiliation operators should have opposite r-dependence, 
one expects that the conjugate pairs are (a+,aft n ), aC n ), and p). Indeed it is 
easy to check that the canonical quantization conditions 

2W 



-5(a - a')r]i 



are satisfied if we impose the following commutation relations: 



X+A 

/x m u mn 1 



A m 5mn j 



(3 J 



A+. 



(4.13) 



(4.14) 



It is then a bit tedious but straightforward to compute the Virasoro operators for the 
£1-^2 system by the expansion 

1 



V L 12 e - in ( T + a ) 

/ j 71 



4a' 



-(^£i + <W + («9 T £ 2 + «9 ff £ 2 ) 



(4.15) 



(One can check that on the right hand side terms which do not have the dependence 
e -m(r+a) a rj cance i.) The result is, with n > 1, 



rl2 
^0 



-PP + E ( a -m a m + a -m a J + T^ 1 ~ V ) 
m>l — 



(4.16) 
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L n = fan ~ fan + H a t a m 

l+m=n;l,m>l 

~\~ ^ t ipt_ m ot m ^ n -\- tt_ m Q? m _|_ n ) , (4.17) 

m>l 

LL 2 n = fa\ -faZ n + a -l a -m 

l+m=n;l,m>l 

+ V (a + a + +a~ a" ). (4.18) 

m>l,m— n>l 

The shift v(l — v)j1 in Lq 2 , which takes exactly the form expected, has been determined 
so that these operators satisfy the usual form of the Virasoro algebra (with central charge 
equal to 2): 



7-12 7-12 

i ^n 



{m-n)L% +n + 5 m+nfi ^{m 3 - m) . (4.19) 



Also note that Lq 2 is hermitian with this shift included. 



Now, since Lq 2 appears to be diagonal in the number operators for the oscillator modes, 



one may expect that the trace over the modes, Trw , readily leads to the expressions in 
( [4. 2| ) . This naive reasoning must however be carefully examined because of the unusual 
hermiticity properties of the oscillators. From the explicit form of £i and £2 given in (|4.9|) 
and ( |4.11| ), one finds that the reality of these fields requires 

("m) 1 " = -aim) ( a n\V = - a -m, ( 4 - 20 ) 

P = P, ? = (4.21) 

This shows that, while (3 and (3 are like a coordinate and its conjugate momentum, 
are interchanged under hermitian conjugation and hence cannot be identified as usual 
oscillators. 

To study the structure of the Hilbert space for these unusual oscillators, let us 
concentrate on a particular level n and consider the system defined by 



[a, a] = X 



b,b 



X* , (4.22) 
a f = -6, S f = -6, (4.23) 
H = aa + bb, (4.24) 

where A is a non-vanishing complex number. The Hamiltonian H is obviously hermitian 
and the system is completely consistent. Suppose we define the vacuum state ^0,0 by 

e^o.o = 0, 0^0,0 = 0. (4.25) 
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Then we can build a general excited state ipm.n by 

</Vn = iV m>n 5 m 6> ,o, (4.26) 

where N m ^ n is a normalization constant. This is clearly an eigenstate of H with the 
eigenvalue mX + nX*, which is complex unless m = n. It is well known that a hermi- 
tian operator can have complex eigenvalues if and only if the norm of the corresponding 



eigenstates vanish |L8|. Indeed we have 

= {bty m , n , VVn,n) = ™A(V>m,n, ^m.n) , (4-27) 

and the norm (VWu VwO vanishes for m ^ n. By a similar manipulation, one can easily 
show that the state which has non- vanishing inner product with i/j m>n is ^ nm . In other 
words, the metric of this Hilbert space is non-diagonaQ If we choose the normalization 
so that (ip n ,m, ipm,n) = 1, then the trace of an operator O should be defined as 

TiO = E (V>n,m, ^Vw) ■ (4.28) 

m>0,n>0 

Hence the partition function is computed as 

TlW H = 0'n,rn,W H tjj rnin ) 
m>0,n>0 



OO 



= e w mX+nX * = E> mA E ^ mA * 

m>0,n>0 \m=0 / \n=0 / 

= (1 - ^Ap^l - w^)- 1 , (4.29) 

which turned out to be identical with the result expected by the naive reasoning. 

To understand the physical meaning of this structure as well as to check the result 
in a different way, it is instructive to form new oscillators which satisfy the conventional 
hermiticity property. For example, if we define A, , B, as 

A = -^=-— (Aa + Xb) (4.30) 



13 This means that the Hilbert space has negative norm states, but this is simply due to the timelike 
nature of £i. 
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where / is the real part of A and is taken positive, they satisfy the standard commutation 



relations 
becomes 



1 



H 



B, 



aa + bb 
|A' 2 < 



1. In terms of these oscillators, the Hamiltonian 



B J B 



2 I A* X 



+; i 



A, 
A* 



A ] B + 



1 



A* 



(4.34) 



Thus, in this representation the metric of the Hilbert space is diagonal but instead the 
Hamiltonian is non-diagonal for complex A. As one can easily check, A and B are, 
respectively, the oscillators which describe the usual particle modes for £i and £2 when 
A becomes real. This means that, except for a special situation where x vanishes, the 
usual particle-like modes are "unstable" and as time goes on they incessantly transform 
into each other. Nevertheless, since the Hamiltonian is hermitian, the total probability is 
conserved and the process is unitary. 

As for the calculation of the trace, although a bit tedious, one can reproduce the 
previous result Q4.29Q using the above non-diagonal Hamiltonian (for example by the use 
of the coherent state method). Finally, applying this to the original problem, we find that 
the contribution of the a^, (m 7^ 0) oscillators to the trace Trw L o 2 is given by 



m+v\—l 1 



W 



m—v\—X 



(4.35) 



m=l 



as anticipated. 

Let us next consider the (3-j3 system, which is characterized by 

■ X 



H 



71 



/3, 



2tt 



(4.36) 
(4.37) 
(4.38) 



Before performing a proper computation of the trace Tru^ = Tre -7rtH ,t = — \nw/n for 
this system, let us describe the essence of the physics by looking at the small x limit. In 
this limit, the boundary conditions for £1 and £2 become almost Neumann and Dirichlet 
respectively. Indeed if we introduce q and p by 



\p,q] 



v 



■ X 

1 

7T 



(4.39) 
(4.40) 
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£j can be written in the form 



6 = ^- + 2a'pr + 0(x 2 ), (4.41) 
X 

6 = -g<x + 0( x 2 ), (4.42) 



with the energy 



H = ^Lqi-a'p 2 . (4.43) 

This corresponds to the following picture: When \ is small, while the a = end of the 
string can fluctuate strictly along the U\ direction, the a = n end may slide along the 
direction which is slightly tilted from u\ into the perpendicular v,2 direction. Thus, the 
string, as a whole without oscillation, can move into the u\ direction almost freely but 
this motion is accompanied by a stretching or shrinking in the 112 direction. This costs 
energy and results in the restoring potential seen in the Hamiltonian. Now from (|4.40|) 
we see that as x tends to vanish, p and q will become independent. Therefore, in this 
limit we get 

Tre -ntH — y - [ dqdpe- ntq2/4a 'e 7Tta ' p2 — , (4.44) 

X J Xlnw 

where the factor in front of the integral compensates the normalization of the commutator 
( |4.40| ) . One can easily check that this agrees with the x ~ > limit of the expression 
1/(1 -w") in §3). 

Let us now compute Tre~ 7TtH for general x- I n order to build a well-defined Hilbert 
space, we shall define a conventional oscillator pair (a, a') by linear combinations of (3 and 
(3. A simple choice that realizes a,a^ = 1, (ay = is 

1 



((3 + ip) , a f = ^_ ((3 - 0) . (4.45) 



(More general ones are unitarily equivalent to this.) Then the Hamiltonian takes the 
form 

H = ^(a 2 -a t2 + l), (4.46) 

27T 

which is not number-diagonal. The relevant trace Tre _7r *^ can nevertheless be computed 
in various ways. One way is to first perform the normal ordering by a technique similar 
to the one explained in the Appendix, which gives 
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e 



-wtH 



5e aa ^e Pa2 e^\ (4.47) 

where 5 = e- ixt/2 (cos X ty 1/2 , (4.48) 

a = ^tan%£, (4.49) 

13 = --sinxtcosxt , (4.50) 

7 = — ln(cosxt). (4-51) 

Then the trace of this expression is easily obtained by using the coherent state method. 
The result is 



Tr e 



-TTtH 



(1 - e 7 ) 2 - 4a/?e 27 l 1/2 (4.52) 



-*X*/2 1 



2i sin ^ 1 — io l 



(4.53) 



Thus we obtain the desired expression. 

We now briefly describe the contributions of the D — 2 transverse components and of 
the ghosts. 

As for £/s, which satisfy the standard Dirichlet boundary conditions, there are no zero 
modes and their contribution therefore is simply f(w)~( D ~ 2 \ 

Finally the ghosts. Since they are not affected by the background, the analysis is 
the same as in the usual case, which was described in detail in the appendix of 
The essential point is the following. Recall that due to the presence on a cylinder of 
a Teichmuller parameter and a conformal symmetry the calculation using the boundary 
states in the closed string channel required an insertion of appropriate zero modes, which 
is equivalent to an insertion of the ghost number operator. When converted into the 
open string channel, the effect of this insertion produces a factor of Inw. In this way, 
one obtains the formula connecting the ghost partition function in the closed channel and 
that in the open channel: 

q^fiq 2 ) 2 = w 1/12 f(w) 2 (-\nw/27t). (4.54) 

This, however, is nothing but the famous modular transformation formula for the (square 
of ) the Dedekind rj function. Since the powers of q and w appearing in this formula 
have already been taken care of by the total intercept of the Virasoro operator, the 
contribution we must add is f(w) 2 (— mto/27r). In particular, this factor of \nw is exactly 
what is needed to cancel the 1/ In w in ( |4.4j ) . 
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Thus, putting everything together, we have reproduced all the factors that appear in 
the scattering amplitude as seen in the open channel. The important lesson we learned 
from this exercise is that non-trivial mixing of particle modes occur in this channel. We 
are intending to understand this phenomenon from the point of view of effective gauge 
theory in a future report. 



5 Discussions 

In this article, we have computed the amplitude for the scattering of two D-particles in 
bosonic string theory, where D-particles themselves are quantized. We employed three 
different methods, namely, the path integral, the boundary state, and the operator for- 
malism in the open string channel, and cross-checked the result. The emphasis was on 
the path integral method, which is conceptually most complete in formulating the prob- 
lem. Especially, it is only with this method that we can clearly grasp the nature of 
the approximation used and compute the corrections in a systematic manner. As far as 
the computation in the lowest order (in the acceleration f{) is concerned, the use of the 
boundary state representation of the interaction vertex appears to be most efficient. On 
the other hand, the operator method in the open string channel reveals the occurence of 
incessant transitions among the excitations of the open string that connect the D-particles 
and perhaps a deeper understanding of this behavior will be important in unravelling the 
connection with the description in terms of the spontaneously broken gauge theory. 

In any event, the fact that the interaction between quantum D-particles through all the 
excitation modes of a string is fully consistent is rather remarkable. This is because in the 
present approach D-particles are no longer backgrounds but are new independent entities 
which can coexist with strings. This raises a puzzling question^]: Is it not true that 
they are supposed to be solitons of string theory and hence they should not represent new 
degrees of freedom ? This may indicate that some important non-perturbative consistency 
condition linking the two is still missing. On the other hand, the conjecture jL9| that the 



D-particles are the Kaluza-Klein modes of 11 dimensional M theory indicates that they 
could indeed be degrees of freedom independent of those of strings. The answer is not yet 
known, but it is certainly an important question in the light of consistency of our result. 

In this article, we have presented the calculation only for bosonic string theory and 
have left out the more important case of superstring. The reason is that in the latter case, 
a satisfactory analysis requires the understanding of not only the scattering of bosonic 



14 



This question was raised by N. Ishibashi during our discussion. 
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D-particles but also that involving fermionic superpartners as well and this appears to 
be non-trivial. In the path integral approach, we need to consider supertrajectories and 
their quantization and this is known to be difficult. This matter is under investigation 
and we hope to be able to report our progress in the near future. 

Nevertheless, if one is content with the amplitude for two bosonic quantum D-particles 
scattering into again two bosonic ones, we can present a well-educated guess. This is due 
to the fact that in the bosonic string case (at least for the lowest order approximation) 
there emerged a simple rule to go from the forward scattering amplitude for infinitely 
heavy D-particles to our general amplitude for quantized D-particles with finite mass 
and it is almost obvious that this rule should continue to hold for the superstring case. 
Assuming this to be the case, all we have to do is to write down the result obtained in 
||, make one simple substitution and a Fourier transformation, and take into account the 
quantum indistinguishability discussed in subsection 2.4. Explicitly, the general on-shell 
amplitude in the open channel representation should read 

•A-super (.Pi i P j ) •A-super^Ply P2i P\i 

p' 2 ) + A super (pi,p 2 ;P2,p' 1 ) , (5.1) 

(5.2) 

where 



A super (PuP2;p[, P ' 2 ) ~ f 'jt-'e-^-^^^ 

° " 1 \2w I 2 ) 



£MglfM0|f)W§). ,5,) 



o=2,3,4 



In this formula, x is as given in (|2.140| ) ~ ( [2.141| ) and e 2 = —e^ = = — 1. We intend to 



discuss in detail how this general formula fits with the super- Yang-Mills description in a 
future communication. 

We finish by making an important observation that the type of amplitude for the 
basic process considered in this paper, no matter how accurately computed in superstring 
theory and no matter how small the string coupling constant is, cannot be relied upon 
if one wishes to know the behavior of the D-particles at very high energy. The difficulty 
stems from the fact that the annulus amplitude does not explicitly depend on the string 
coupling constant g s . Thus, one can insert "annulus interactions" anywhere one likes and 
produce infinitely many diagrams, such as the ones dipicted in Fig.4, which contribute at 
the same order in g s . 
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Fig. 4 Example of diagrams at the same order in the string coupling as the basic one. 

At energy much smaller than the D-particle mass, the diagrams with extra corrections are 
expected to be suppressed by powers of E 2 /m 2 due to the presence of internal propagators 
of D-particles. But as the energy becomes comparable to m all these diagrams would 
contribute equally to the amplitude and computation appears to become intractable. It 
would however be extremely interesting if one can devise a method to extract the leading 
high energy behavior out of these diagrams just as the similar study for ordinary strings 
20fl stimulated a lot of thinking about the degrees of freedom of string theory at short 



distance. 
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Appendix: Normal Ordering Formula 

Let (tZj, a\) and (a*, a\) be two independet sets of oscillators satisfying the usual (anti- 

= 5ij. Throughout, upper (lower) sign refers 



) commutation relations 



(2^, dj 



to the anti-commuting (commuting) case. Consider the bilinears 



a 



aAa , (3 = a'Ba' 



(A.l) 



where aAa = J2ij cuAijdj, etc.. By extending the idea of frTfl , we will derive a formula 
which gives the fully normal ordered form of e a e 13 . 
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Let A be a parameter and we seek the result in the following form: 

x(X) = o t X(A)a t , (A.3) 
y (X) = a Y(X)a, (A.4) 
z(X) = rfZ(X)a + a^W(X)a. (A.5) 

Here X(X), Y(X), Z(X), W(X) are matrices and 5(A) is a prefactor, to be determined. Let 
us define 

Xi(A) = 6(X)e*We^\ (A.6) 
X2 (A) = e Ao e A V* (A) . (A.7) 

Then what we want to obtain is Xi{X) = X2{X) and this is achieved uniquely if we can 
satisfy 

(i) Xi(0) = X 2 (0), (A.8) 
(«) xr 1 (A)x'i(A) = X 2 - 1 (A) X / 2 (A), (A.9) 

where x'\ = d-Xi/dX etc.. Writing out the condition (it) explicitly, we get 

j + Xl 1 x'(X)xi + y'(X) = X2 1 aX2 + X2 l e Xa Pe- Xa X 2-z'(X), (A.10) 

where we have made an assumption [z(X), z'(X)] = 0. Its validity will be aposteriori justi- 
fied. Evaluation of the quantities appearing in ( |A.10| ) is a bit tedious but straightforward. 
Then, by equating the coefficients of the same operator structure, we get the following 5 
differential equations: 

(a) ^ = XTt(BA) + Tr(AT) , (A.ll) 
o 

(b) X' = e z (l t X 2 BA) Be wT , (A. 12) 

(c) Y' =f YX'Y = e~ wT Ae~ z , (A.13) 

(d) X'F = tZ' - Xe z BAe- z , (A. 14) 

(e) YX' = W' T - Xe~ wT ABe wT . (A.15) 

From the structure of these equations, it is consistent to assume that X'Y, Z, X'A, BY 
are functions only of BA, while YX' and W T depend only on AB. Let us take the trace 
of (d) and (e). We get 

Tr(XT) = TrZ' — ATr(£M) 

= TtW — XTi(BA) . (A. 16) 
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Thus we must have TiZ' = TrW. Substituting ( |A.16 ) into (a) and integrating, we get 



S = e TtZ = e TrW . (A.17) 

This result, together with the previously mentioned dependence on AB and BA, prompts 
us to make a postulate 

Z = f(BA) , W T = f(AB), (A.18) 

where f(x) is some power-expandable function. Then since Bf(AB) = f(BA)B holds, 
we have 

BW T = ZB, W T A = AZ. (A. 19) 

Using these postulates, the equations (b) ~ (d) simplify to 

(b) X' = e 2Z (l =f A 2 £) B , (A.20) 



(c) Y' t YX'Y = Ae~ 2Z , (A.21) 

(d) X'Y = tZ' - A£ , (A.22) 
where ^ = BA . (A.23) 

To solve them, first multiply (b) from right by Y. We get 

X'Y = e 2Z (l =F A 2 ^) 7] , (A.24) 
where r] = BY. (A.25) 

Equating this with the RHS of (d), we can solve for rj in terms of £: 

V = (ItA^-V^t^-AC). (A.26) 

Next multiply (c) from left by B and substitute (d) for X'Y. One then obtains 

v ' Tv ( T Z'-\0 = ie' 2Z . (A.27) 

We may now eliminate r\ by substituting ( |A.26| ). After some calculation we get a differ- 
ential equation for Z of the form 

(\^\ 2 i)(Z' 2 - Z") = ±2A£Z'±2£. (A.28) 

The solution of this equation satisfying the proper boundary condition Z(0) = is given 
by 

Z = -ln(l±A 2 0- (A.29) 
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We then get from ( A. . 1 7| ) and the postulate ( A.l 



1. 



±1/2 



5 = exp^±-Trln(l±A 2 0J = [det(l ± A^)J > ( A - 30 ) 

W = f{B T A T ) = - In (l ± X 2 B T A T ) . (A.31) 

The matrix X, satisfying X(0) = 0, is obtained by integrating (b). We get 

x = [ ' du JTT^ B = ttxh b ■ (A - 32 > 



To get Y, we recall Y — B 1 rj, where rj is given by (|A.26|) . One immediately gets 



Y = AA(1±A 2 0- (A.33) 

The remaining equation (e), which now reads 

YX' = tW' T -XAB, (AM) 

is automatically satisfied. Also one can easily check the validity of the assumption made 
earlier, namely [z(X), z'(X)] = 0. Since the solution of the system of equations with 
appropriate boundary conditions is unique, we can justify all the postulates made above. 
Putting all together, we get the normal-ordering formula quoted in the text: 

= [det(l±BA)] ±1 e a ^ 1±BA ^ lBa ' 

-e aA(l±BA)a e -aUn{l±BA)a e -aHn{l±B T A T )a gc^ 

For completeness, we exhibit the formula for the case of one set of oscillators, which 
can be obtained in a similar manner: 

giaAdgiat Bat 

= [det(l± J BA)] ±1/2 e^ t ( 1±B ^ lBat 

. e \aA{l±BA)a e -aUn(l±BA)a _ (A.36) 

In this case, the matrix A and B are symmetric (anti-symmetric) if the oscillators are 
commuting (anti-commuting). 
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